Residual Resistance in a 2DES: A Phenomenological Approach. 
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We consider a simple phenomenological model of a semiconductor with absolute negative con- 
ductance in a magnetic field. We find the form of the domains of the electric field and current which 
arise as a result of an instability of a uniform state. We show that in both Corbino disc and Hall bar 
samples the residual conductance and resistance are negative and exponentially small; they decrease 
exponentially with increasing length L Xty . 



I. INTRODUCTION 



An interesting effect recently observed is that the resistance of a microwave irradiated two dimensional system 
(2DES) in a magnetic field B drops almost to zero in certain intervals of B [1,3]. In magnetic fields too low for 
the observation of Shubnikov-de Haas oscillations (B < 0.1T) the resistance R xx of a 2DES with mobility ~ 3 x 
10 e V/ cm. sec 2 begins to oscillate in the presence of microwave irradiation [2]. In samples with higher mobility 
(~ 2 x 10' V/ cm. sec 2 ) the resistance oscillations become more pronounced and the resistance drops almost to zero in 
CO . some intervals of B [1,3] lower than the fields at which Shubnikov-de Haas oscillations begin to appear. Sometimes 
the states of the 2DES with almost zero resistance are called the zero resistance state (ZRS). These states have 
been observed in experiments on Hall bars [3,1]. Similar measurements have been recently carried out on Corbino 
■ discs [4]. In this case a zero conductance state rather than the ZRS has been observed in the same intervals of the 
, applied magnetic field B. These experiments have drawn considerable attention, and during a short period of time 
many papers appeared where possible mechanisms for the observed effects were suggested [5-16]. The most adequate 
mechanism seems to be the one based on the idea of absolute negative conductance (ANC) [6-10,12,13,15,16]. It was 
shown that a microwave irradiated 2DES in a strong magnetic field (ui c t > 1, where u> c is the cyclotron frequency, r is 
the momentum relaxation time) may have absolute negative conductance a xx . In some papers the calculations were 
carried out for a small bias electric field E (a linear response). In Refs. [7,16] the dependence a xx (E) was discussed and 
the conclusion was made that the dependence of a xx {E)E x on E x may have a N-shape form (see Fig.l). The authors 
of Ref. [8] assumed that the function E x = p xx (j x )j x has the N-shape form and showed that the states with negative 
differential conductance (NDC) are unstable with respect to small perturbations. As a result of this instability, a 
nonuniform distribution of the current density is assumed to arise and current domains (or current filaments) should 
O ■ appear in the sample. In this state the field E x is zero and a change in the total current leads for example to an 
enlargement of the domain with positive current and to a shrinking of the domain with negative current (the ZRS). 

The possibility of a microwave induced ANC in semiconductors in a quantizing magnetic field was suggested and 
discussed some decades ago [24-26] . Around the same time phenomena in semiconductors with N- and S-shape I(V) 
characteristic were studied inensively both experimentally and theoretically (see the review article [18] and the book 
[19] and references therein). In particiular, it was shown that in semiconductors with the N-shape /(V) characteristic 
generally speaking moving domains of the electric field arise as a result of the instability of a state with the NDC. In 
semiconductors with the S-shape /(V) characteristics, domains (filaments) of the current density arise as a result of 
the instability. In the presence of the electric or current domains the I(V) characteristic changes drastically. In the 
first case an almost horizontal section appears on the I(V) characteristics (zero differential conductance), whereas an 
almost vertical section appears on the I(V) curve in the second case (zero differential resistance). If a system with 
S- or N-shape I(V) characteristic is placed in a magnetic field, the situation becomes more complicated: the form of 
the I(V) curve now depends on the measurement geometry (Hall bar or Corbino disc) [12] and both electric field and 
current domains may coexist. For example, if one assumes the N-shape form for j x = o~ xx (E x )E x dependence on E x 
in a Corbino disc (E y = 0), then in the Hall bar measurements with increasing magnetic field the dependence j x (E x ) 
acquires a complicated form becoming S-shaped at ui c t > 1 [12]. 

Strictly speaking, equations for macroscopic parameters (electric field, electron density etc) corresponding to ex- 
perimental samples [1,3,4] are not yet derived. In the recent paper [16] a kinetic equation for the distribution function 
is derived which, in principle, could be applied to describe nonhomogeneous states. However this task is very diffi- 
cult. All previous publications on these nonhomogeneous states are based on equations for macroscopic parameters. 
Therefore, although the dependence of the ZRS or ZCS on the applied magnetic field B obtained theoretically [6,7] 
qualitatively agrees with experimental results [1,3,4], it is difficult to convincingly conclude that the observed ZRS or 
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ZCS are the result of instability of an uniform state in the 2DES with the S- or N-shape I(V) curve. Nevertheless 
even a phenomenological model of the S- or N-shape characteristics allows one to make certain conclusions about 
the observed I(V) curves. In this paper we admit a phenomenological model assuming a N-shape I(V) curve in 
the Corbino disc geometry and calculate the residual resistance (Hall bar experiments) or conductance (Corbino disc 
experiments). We show that in our model this resistance (conductance) depends exponentially on the ratio of the 
width of the sample to the width of the domain wall. Although our model does not correspond to the parameters 
of the real samples used in experiments [1,3,4], it has a physical meaning and allows us to make general statements 
about the stability of the system, the form of electric and current domains and the residual resistance of the system in 
a nonuniform state. In particular using model one can answer the question whether the state observed in experiments 
is dissipationless (as it happens in superconductivity and was anticipated in Ref. [1]) or not. 



II. MODEL 

As in Ref. [12], we assume that a xx depends on the electric field E — \J E^ + E^ so that in the absence of E y 

(the Corbino disc) the dependence j x — a xx (E)E x on E x has an N-shape form (see Fig.l). Such a dependence was 
discussed in Ref. [7] (see also [16]). For example, this dependence may be obtained if we take for a xx (E) the expression 

axx{E)=a °°-i+(*oM(E/E r (1) 



where (Too is the conductance in strong fields (E >> E ), a = —a xx (0), and a xx (E) equals zero at the field E - 
The electric field Eq as well as the conductivity <jq should be determined from a microscopic theory (see 
Refs. [7,16] and [17]). For example, according to Ref. [17] on the order of magnitude Eq > 20 V/cm. In 
order to describe the form of domains, one has to know a microscopic mechanism which determines the form of the 
I{V) curve. For example in Refs. [18,23] a superheating mechanism of the S-shape current-voltage characteristic was 
considered. In this case the width of the current domains is determined by the energy relaxation length. We assume 
that in our model a characteristic length is determined by the screening length (such a model was used to describe 
the Gunn effect [18,19]). Therefore the current density can be written as 

j = <tE - et)Vn + (e/4ir)dE/dt (2) 



where a is the conductivity tensor with diagonal components equal to a xx = a yy = a and off-diagonal ones equal 
to a xy = —a yx = an , D is the diffusion coefficient tensor which has a similar structure. The last term is the 
displacement current. We assume that only a depends on the field E. The components an and D, Dh are assumed 
to be independent of E. At first glance the assumption about independence of the diffusion coefficient 
on E violates the Einstein relation between the mobility and diffusion constant. However this relation 
holds only for systems in equilibrium. This is not the case for the system under consideration. Of 
course the diffusion coefficient D depends on E in some (unknown in our phenomenological appoach) 
way. This dependence can be taken into account if we introduce a new "potential" in Eq.(6) or a 
new electric field. Therefore a corresponding analysis for the case of the energy dependent diffusion 
coefficient D(E) can be carried out in a similar way as for the constant D. For simlicity we ignore this 
dependence because it does not lead to qualitatively new results. Note that an analogous problem was 
discussed in the theory of the Gunn effect (see the reviews [18,19] and references therein). 

The electron concentration n is related to the electric field via the Poisson equation (instead of -e we write +c 
making a proper choice of the electric field direction) 

VE= (47re/e)(n-n ) (3) 

We also assume that the thickness of the 2DES d is larger (or of the order) than the screening length l scr . In this 
case in order to find a relation between n and E, one has to solve Eq.(3) inside the sample. Otherwise one has to 
solve an equation for the electric field (or potential) outside the sample taking into account corresponding boundary 
conditions. In the case l scr < d Eqs.(2,3) yield correct results at least qualitatively. The main drawback of our model 
is the assumption of a local relationship between the current density j and electric field E which is true if the mean 
free path I is shorter than a characteristic length of the problem (in our case the screening length) . In the experiments 
one has the opposite relationship between l scr and I. 
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FIG. 1. The assumed I(V) characteristic in the Corbino disc measurements (solid line) and the corresponding conductivity 
a xx as a function of ^(dashed line). Here a m i n = <ro and a max — a^. 

III. THE CORBINO DISC 

Let us consider first the simplest case of the Corbino disc geometry when E y = (see Fig. 2a). 
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FIG. 2. The Corbino disc (a) and Hall bar (b) samples (schematically). The region with diagonal stripes and the cross-hatch 
region mean domains with opposite directions of the current 



Eqs.(2,3) can be rewritten in this case as follows 

jx = vE x - eDd x n + (e/4n)8E x /dt 



(4) 



d x E = (47re/e)(n-n ) (5) 

These equations practically coincide with equations used in the theory of the Gunn effect [18,19]. Linearizing these 
equations with respect to small fluctuations E(x, t) — E x q + 8E x {x, t) * exp(icj£ — ikx), one can easily show that states 
with negative differential conductance are unstable: iu> = —djo(E x )/dE x — k 2 D, where jo = cr xx (E x )E x is the /(V) 
characteristic in an uniform case. This dispersion relation shows that the states with negative differential conductance 
are unstable with respect to long- wave perturbations. The characteristic length of the problem (the width of a domain 
wall) is determined by the relation: l c h = y/D/ fa^J, where ad = 8jo(E x )/dE x is the differential conductance. The 
conclusion on the instability of the states with a negative differential conductance was made a long time ago (see Refs 
[20-22]). As a result of the instability, domains of the electric field (in the x-direction) arise in the sample. In the 
y-direction these domains may be regarded as the current density filaments: different Hall currents flow in domains 
with opposite directions of the electric fields. In order to find the form of the domains of the electric field E, we 
exclude n from Eqs.(2,3) and arrive at the equation 

3 2 E N 8U dE N . dE N 

-JT2~ + 7^ + 3no{E n )- =0 (6) 
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Here En — E x /Eq is the normalised electric field, xn — x/l scr , l scr — J De/AniJoc is the characteristic screening 
length, 7 = eE Q /4:Tren a l scr , jjvo(-Ejv) = ja(E N )/ Eoa^, t N = t/t ch , t ch = en {) l scr j o^Eq. The "potential" U is defined 
as 

/•En 

U(E N ) = / dE Nx \j Nb - 3no(E Nx )] (7) 
Jo 

where jNb = jb/^ocEo is the normalised bias current, Jno(En) = E a{E N ) / a^E^ is the normalised I(V) curve 
(see Fig.l). The electron density n in an uniform case is the 3-dimensional electron concentration: it is related to 
the 2-dimensional density ri2o via no = 11,20/ d, where d is the thickness of the 2DES. We neglect the last term in 
Eq.(6) assuming that 7 is small. Let us estimate 7 for the samples used for example in the experiments [1,3,4]: for 
^20/^2 ~ 3 10 11 cmT 2 /lscr, we get 7 = Eq/E c h, where E c h ~ 4 10 4 V/cm. Obviously the field Eo is much smaller 
than the very large value of E c h- To find the form of stationary domains of the electric field, we neglect the last term 
in Eq.(6). In the following the analysis is similar to that presented in Rcf. [18]. We integrate Eq.(6) and obtain the 
"energy" conservation law which describes the phase trajectories 

(1/2)E* + U{E N ) = U(E N1 ) (8) 

The integration constant U(Eni) determines the phase trajectory (see Fig. 3). Its choice depends on concrete 
boundary conditions. We assume that at the boundaries x = ±L the spatial derivates E' N = OEn/Oxn is zero (the 
final result does not depend qualitatively on particular boundary conditions). We are interested in a solution of the 
form of two domains with opposite electric fields. The field E N \ is connected with the length of the sample L x by 
the relation 

/Eni 
dE N /E' N (9) 
-E N2 



where = L x /l scr , L x is the distance between the outer and inner radious of the Corbino disc. The constant E^2 
is connected with E^i by the equation 

U(-E N2 ) = U(E N1 ) (10) 
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FIG. 3. The phase trajectories. The solid line shows a trajectory which describes the electric field distribution in the 
considered finite sample in the presence of the bias current jbN (two domains with opposite directions of the fields). The 
maximum fields in these domains are equal to —En2 and Eni- The dash-dotted curve shows a trajectory which corresponds to 
an infinite domain with a positive field and a finite domain with a negative field (soliton). Two infinte domains with opposite 
electric fields in the absence of the bias current are described by the dashed trajectory. 

Our goal is to find a relation between the averaged field in the sample E^ with the bias current jb- The relation 
jb(Eb) is the form of the current-voltage characteristic in the presence of domains. The averaged field E b = EnE is 
given by 

dEx En/E' n (11) 

-En2 
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The equations (9,10,11) determines the jb(Eb) characteristic, or the residual conductance in the presence of domains. 
The function E n (En) in these equations is shown in Fig. 3 (solid line). It is clear that if the characteristic width of 
the domain wall (in our model l scr ) is much less than the length L x , the trajectory E n (En) should be close to the 
dash-dotted trajectory (the soliton trajectory). We calculate the residual conductance G res assuming that l scr « L x . 
In this case the fields ±Eni,2 inside domains are close to ±1, the normalised bias current jbN is very small and in the 
main approximation the integrals in equations (9,10,11) can be calculated analytically. Taking this into account, we 
find from Eqs. (9, 10, 11) 

5E% 2 w 5El + 46 (12) 
(5E N1 -b)(SE N2 + b) = cxp(-V2^L N ) (13) 

Bw = (1/i „^ )ln (g=±| (14) 

Here SE N1 = 1 - E Nl , SE N2 = -E N2 + 1, b = -jwv/2a, a = (l/2)dj N {E N )/dE N \ EN =i = 
(l/2)5<7jv(.Ejv)/9-EjvJ.Ej V =i- Note that in our model the average field Em is positive if the bias current jbN is nega- 
tive, i.e. the resudial conductance is negative. We calculate the residual conductance in the limit of small currents 
jbN « 1- If the condition 



is satisfied, we get from Eqs. (12, 13, 14) 



SE 2 N1 » b (15) 

SE N 2 = SE N i [1 + 2b/SE N1 ] (16) 
SE m = exp(-V2aXAr/2) (17) 

E N = 2bcxp(V2^L N )/(V2aL N ) (18) 
Therefore the residual conductance has the form (we restore the dimensional units) 

G res = -aoo(aV2aL x /l scr ) exp(—V2aL x /l scr ) (19) 

We see that b as well as 5Eni,2 are exponentially small. Therefore the assumed smallncss of 5En2 is justified. The 
formula (19) is valid if the condition (15) is fulfilled. This condition can be represented in the form 

E « E Q {l scr /L)V2a (20) 

The formula (19) shows that the residual conductance is small (we assumed that L x » l scr ) and exponentially 
decreases with increasing length L x . Note that in a model more appropriate to the real experimcmcnts [3,1] the 
characteristic length might differ from the screening length l scr . If E exceeds the value of the characteristic field in 
the rhs of Eq.(20), the residual conductance increases exponentially. Note that the obtained nonhomogeneous state 
with a negative differential conductance is stable if the total voltage is fixed (see Ref. [18]). 

IV. THE HALL BAR 

In this Section we consider the case of the Hall bar geometry (see Fig. 2b). In this case we have to write the equations 
for the currents j XyV having in mind that all quantities depend on y 

j x = oE x + a H E y - eD H d y n + {e/4n)dE x /8t (21) 
j y = oE y - <j H E x - eDdy-n + (e/A^dEy/dt (22) 
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Taking into account the Poisson equation (3) and the fact that j y = 0, we obtain 

j b W G R E y - {eD H /A-K)dyyEy 



(23) 



h « (7(£7 W )£7 B - {eD/4*)dyyE v (24) 

Here = j x (y) is the bias current density which depends on y, j\ = grE x . The conductivity g(E) w cr(E y ) because 
we assume that 07/ >> a, and therefore £^ >> E x . As before we neglect small terms of the order 7. The form of 
domains is determined by Eq.(24) that coincides with Eq.(6) if the latter small term in Eq.(6) is neglected. We note 
the electric field component E x does not depend on y (this follows from the Maxwell equation in the stationary state 
which we are interested in: V x E — 0). The current jb is y-dependent; it has opposite directions in different electric 
domains. Our aim is to find a relation between the average current j b and the electric field E x , or the inverse relation 
E x Qb). From Eq.(23) we find 

lb = °H~Ey (25) 

Calculating a relation between j\ and E y from Eq.(24), as we did in the previous Section, we get 

ji = -G res E y (26) 
Combining Eq.(25) and Eq.(26), we finally find 

E x = -(G res /a H )] b (27) 

This result means that the residual resistance R res is equal to R res = {G res / g h ), i.e., the residual resistance also 
is exponentially small in the Hall bar samples measurements. Note that in the expression for G re s (19) the length 
L x should be replaced by L y . The I(V) characteristics for a 2DES with two current (electric field) domains is shown 
in Fig.3 of Ref. [12]. 



V. CONCLUSIONS 



Using a simple phenomenological model, we have calculated the residual conductance (resistance) of a 2DES. We 
have shown that in measurements on both the Corbino disc and the Hall bar the residual conductance and resistance 
are exponentially small; they decrease exponentially with increasing length L x ^ y . In our model the residual conductance 
and resistance are negative. A direct comparision of our results to experimental data meets difficulties because, as 
we mentioned before, parameters of the model do not correspond to real samples. Nevertheless one can note that the 
observation of a small negative resistance was observed in Ref. [27]. It would be of interest to measure the dependence 
of R res (in Hall bars) or G res (in Corbino discs) on the length L y or L x . This can shed light on the applicability of a 
model based on the negative conductance to the phenomenon observed in the experiments [1,3,4]. 

Note that we considered a limiting case gr » a. For an arbitrary ratio gr/g a linear analysis of the 
superheating instability in semiconductors with the S- or N shape I(V) characteristics has been carried 
out by Kogan [30]. It was shown that, as it happens in plasma [31], the superheating instability results 
in the appearance of oblique current filaments (electric field domains). Only in the limit gr/g » 1 
these filaments are parallel to the bias current. Analysis of oblique domains is much more difficult 
problem as one needs to solve nonlinear equations in two dimensions with corresponding boundary 
conditions. It is clear however that if the voltage difference V is measured between peripheral contacts, 
the Hall voltage Vr also contributes to V in the case of oblique domains. Since Vr changes sign with 
changing the magnetic field direction, the voltage V also can change sign. Therefore the contribution 
of the Hall voltage Vr to V might be a reason of the sign change of V with changing the magnetic field 
direction observed in recent experiments [27]. 

Note added: After the preparation of this manuscript [28] we became aware of Ref. [29] in which similar ideas 
were elaborated. 
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